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Massive gravity in the weak field limit is described by the Fierz-Pauli theory with 5 degrees 
of freedom in four dimensions. In this theory, we calculate the gravitomagnetic effects (potential 
energy) between two point-like, spinning sources that also orbit around each other in the limit 
where the spins and the velocities are small. Spin-spin, spin-orbit and orbit-orbit interactions 
in massive gravity theory have rather remarkable, discrete differences from their counterparts in 
General Relativity. Our computation is applicable for large distances, for example, for interaction 
between galaxies or galaxy clusters where massive gravity is expected to play a role. We also 
extend the computations to quadratic gravity theories in four dimensions and find the lowest order 
gravitomagnetic effects and show that at small separations quadratic gravity behaves differently 
than General Relativity. 


I. INTRODUCTION 


One of the simplest modifications of General Relativity (GR) with improved infrared (IR) behavior that could 
possibly explain the observed accelerated expansion of the Universe is Massive Gravity (mGR), a topic which has 
been around since the work of Fierz and Pauli in 1939 [T] but received a revived interest in the past decade for various 
reasons and in various forms (See the excellent reviews on the subject |11[5])- At the lowest order, in four dimensions, 
the Lagrangian density is that of the Fierz-Pauli (FP) theory 
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( 1 ) 


which propagates a non-ghost, non-tachyonic massive (with mass rrig) spin-2 particle with all 5 helicity modes in flat 
backgrounds. By construction, at the level of the Lagrangian, or at the level of the field equations, as rUg -T 0, GR is 
smoothly recovered. But it is well-known that once the Newtonian potential is computed between two static sources 
or when deflection of light is computed as it passes a static source, one of these two results can be matched to the 
GR value by redefining the Newton’s constant but the other one does not smoothly reproduce the Newtonian (or the 
GR) limit. For example the potential energy in mGR between two static sources reads 


UmGR 
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( 2 ) 


and as mg -T 0, one has the famous van Dam-Veltman-Zakharov (vdVZ) [HIS] discontinuity which has been known 
since early 1970s. All this is well established in the literature. But, what is rather interesting is that, only recently 
|S|, another discontinuity between GR and mGR was found in the spin orientations of spinning point masses (such 
as two widely separated stars, galaxies or galaxy clusters). In GR, the spin-spin interaction (when other interactions 
such as tidal forces etc are neglected) can be computed as 


j ^spin-spin 
^GR 



Jl . J2 


3Ji . r J 2 . f 


(3) 


where Ji and J 2 are the spins of the sources and r is the radial vector between them. (Here, spin refers to rotation 
about an axis passing through the object which we approximately take as point-like.). On the other hand, in mGR, 
the same quantity can be found as |5] 


^ ^spin-spin 
^mGR 


Ge (1 -I- mgr + m'^r'^) 


Jl . J 2 — 3 Jl . f J 2 . r 


(1 -I- mgr + ^rngT^) 

1 + mgV + m^r"^ 


(4) 
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Figure 1: Minimum energy configuration in GR, as long as weak field limit is applicable: spins are anti-parallel to each other, 
so the spin is minimized. 


Observe that as rUg —>■ 0, exponential decay part is expected since massive gravity is weaker at large distances and 



Figure 2: Minimum energy configuration in massive gravity: for distances r > 1.62/mg, spins are parallel and point in the 
direction perpendicular to the axis joining the sources. For small separations mGR reduces to GR. 

even perhaps, one could expect, just like the Newtonian potential case, a discrete numerical difference in the overall 
factor of the potential. But, what is highly surprising is the discrete numerical difference between the two spin-spin 
terms in Q: As r —>■ 00 , one has 
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Ji . J 2 - Ji . r J 2 . f 
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which differs discretely from that of GR (§. But that is not the whole surprise : As we expect from the Newtonian 
attraction such a discrete discontinuity might arise. The real surprise is that this numerical difference leads to 
different observable spin orientations in these two theories. While GR favors anti-parallel spins ( as shown in Figure 
1), minimizing the total spin of the system; massive GR favors parallel spins perpendicular to the line joining the 
sources, maximizing the total spin of the system ( as shown in Figure 2) . The fact that a tiny mass leads to such a 
remarkable effect is rather amazing. A detailed analysis in [5], actually shows that for a given graviton mass rUg, for 

1 

distances that satisfy nig r > —— « 1.62 the spin of the system is maximized as shown in Figure 2. On the other 


hand, for distances that satisfy nig r < 


1-k V5 


, massive gravity agrees with GR and spin of the system is minimized. 
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At the point rUg r* 1.62, the relative coefficient between the two terms in Q becomes —2. Therefore, as expected 
at the large separations in the Universe, mGR can show its effects on the orientations of galaxies or galaxy clusters. 
In the Gonclusion section, we will point out to some possible observations along these lines. 

Encouraged by these interesting results, in this work, we continue our analysis on the other and higher order 
gravitomagnetic effects in massive gravity as well as higher derivative gravity theories. We shall assume that the two 
spinning sources also have velocities, hence orbital motion which will yield additional effects that could probably lead 
to interesting differences between GR and mGR. 

Our computations in quadratic gravity are in some sense academic in nature since the differences will arise in UV 
not in IR: At large separations, clearly one expects GR to dominate over the quadratic terms, but at small separations, 
quadratic terms are more effective. 

Let us say a few words about how we shall proceed to compute the promised effects. The "canonical" way of 
finding gravitomagnetic effects in a given gravity theory at the weak field, small velocity, small spin limit is to find the 
linearized Kerr solution of the theory and to compute the "motion" of a spinning, orbiting test particle in this geometry. 
But this is a rather long and cumbersome procedure, instead of that we shall employ a more clean-cut method, that is, 
we will compute the tree-level scattering amplitude due to one graviton exchange between two covariantly conserved 
sources. Details of this procedure has been laid out in [ 7 ] but it pays to recapitulate here for the sake of completeness. 

The lay-out of the paper is as follows: In section II, and III we show how the gravitomagnetic potential energy 
including all the terms up to desired order can be computed from the tree-level scattering amplitude. In section IV, 
which is the bulk of the paper, we give a detailed derivation of the potential energy for the moving, spinning masses 
in both GR and mGR. Some of the details of the computations are relegated to the Appendix. In section V, we 
calculate the analogous expressions for quadratic gravity (with no explicit Fierz-Pauli mass term). 


II. POTENTIAL ENERGY FROM GRAVITON EXGHANGE:TREE-LEVEL SGATTERING 

To relate the potential energy at the desired order to the scattering amplitude, let us follow [5] and compute the 
vacuum to vacuum transition amplitude in the path integral formalism between two sources: 


(0|e-*^*|0) 



( 6 ) 


where t is a large time that will drop at the end and S[h,T] is the linearized action about a background {gg,i/) in a 
generic gravity model which reads 

S[h,T] = j d^x^g[ - (7) 

From this follow the linearized field equations 




( 8 ) 


Govariant conservation of leads to = 0. Let us employ the background field method to recast the 

path-integral in a more manageable form. For this purpose, suppose hg,i, satisfies (j^, then make a change of variables 
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in the path-integral as —>■ + h^i, which does not change the measure, but shifts the action to the decoupled 

form 

S[h,T] = I d^x^{ - (9) 

The second term is free of so we can move it out of the path-integral and the first term, being independent of 
simply rescales the normalization factor, yielding 


W^[r] (10) 

Equation ([^ is of the form 

= |T“^(:r), (11) 

where is a self-adjoint operator whose Green’s function is defined as 


^fiuotfuG ^ ap (^5 ^ ) — 2 T 9ppQj^f^ 5{^x x ). 

Therefore the particular solution of can be formally written as 


( 12 ) 


hpu{x)='^ J d^x'y^G^„aij{x,x')T°‘^{x'). (13) 

Therefore one has the usual source-source interaction 

1E[T] = AfeT f d°x'y/-g{x)yy-g(x')T'^‘'(x)G^^^p{x,x')T°^[x') ^ 

from which we can read the potential energy (up to an irrelevant constant) as 

U = -^ J d^xd^x'y/-g{x)y/-g{x')T^'' (x) (x,x') T°‘^ (x'). (15) 

We have kept the discussion to the H-dimensional, for _D = 3 -I- 1, we have k = 16ttG. Figure 3 represents the 
interaction. 


III. COMPUTATION OF THE POTENTIAL ENERGY 


We shall compute the mentioned gravitomagnetic effects in three different theories: General Relativity, (linearized) 
massive gravity and quadratic gravity. To be able to employ (151 we need the propagators of these theories. All 
these theories have different propagators which require a lot of work to carry out the full computation. To reduce the 
amount of computation, we shall use the "master" scattering amplitude that was found in [7], which upon taking the 
corresponding limits, generates the relevant scattering amplitudes in these three theories. There is a caveat though, 
to be able to smoothly reproduce the GR’s and mGR’s scattering amplitudes from one single expression, one must 
introduce o provisional cosmological constant which is set to zero before rrig —> 0 0,110]. The relevant action is. 


S = 


+ J d^X I “ ^ + ^matter | , 



( 16 ) 
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Figure 3: Tree level scattering digram, between two conserved sources. From this diagram, we compute the potential energy. 


and the tree-level one graviton exchange in this theory is 0 

(xO |(/3Q + a)(Af - ^ 


AUt = 2 


D- 1 


r {{pD + a){D + 


4A , m-i 
D-2^ ~ 


T 


(x) 


(17) 


4A 


{D-2){D-IY 
2 


, I - - 4A TTln 

T' (/3n + a)C+— )-^ 


{D-2){D-l) 





which looks rather cumbersome, but bear in mind that it encompasses all the information one needs for the theories 
that we are interested in. □ is the d’Alembertian in the background. The parameters that appear in this scattering 
amplitude are defined as 

1 4AD 4A „ 4A (L» - 3) (D - 4) 


K D-2 D-l 
A{D -l)a + DP 


[D-l) {D-2) 


0 = ^ (18) 

where the effective cosmological constants can be found from -|- /A^ = 0. For GR, in flat space the proper limits 

rrii , 

are 0, a = /3 = 7 = 0, a= - which yield 


K 


2Ut = + —7_T'(52)-ir, 


D-2 

where we have suppressed the integral signs, but more explicitly one has 

n 


2Ut = -K J d°x J d^x’T^^ {x')G{x,x')Tf^'' J j d^x'T {x') G {x,x')T (x) 


(19) 


( 20 ) 
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where we have employed the scalar Green’s function which is easier to handle 

d^G {x,x') =—5 {x,x'), (21) 

with the flat space d’Alembertian 9^ = —d^ + V^. In the explicit computations below, we use the retarded Green’s 
function which are easily found as 

= Gr {x, x') = J t ')], (22) 

47r 2 


for the massless case and 


_ fyiiL) — 

(a^ - ml) = Gr {x, x') = ^ Kd^ (r mg) S[r- {t- t')] , (23) 

(27r)~ " 

for the massive case. Here K^, is the modified Bessel function of the second kind. Up to this point we have kept 
the discussion to be in generic D dimensions, but from now on we shall stick to the D — 3 + 1 dimensional case. 
In [S], spin-spin interaction was found in D dimensions but because the spin is a higher rank tensor in more than 4 
dimensions, the expressions become quite non-trivial, even if the computations are straightforward. 

Before we close this section, let us digress a little bit and compute the interaction of two massless sources (say 
photons) moving at the speed of light, the ultra relativistic limit that we shall avoid in the rest of this work, and 
reproduce the old, amusing result of m, that two anti-parallel moving light beams attract each other, while two 
parallel moving beams do not interact at all (note that [S] has a nice discussion on this issue). Let the four velocities 
be 


< = (1,0,0,1), w^ = (l,0,0,a), (24) 

where cr = 1 refers to parallel motion and tr = — I anti-parallel motion. The energy-momentum tensor for each photon 
is 


= Eu'^u'', 


(25) 


with T = 0. Vanishing of the trace immediately says that the result will be the same in GR and mGR. Then one has, 

4Ut = -2KT'^^ 

= -2kT;,o - 2kTL TG (26) 

Upon carrying out the integrals one arrives at 






t, 


which yield the potential energy between these two photons as 


U = 


-2GE, Eo 




(27) 


(28) 


So clearly photons that move parallel to each other do not see each other {U = 0) at this level of approximation 
(namely the weak held limit and neglecting the spins). As noted this is also valid in mGR . Gan this amusing result 
have observable consequences? Imagine two photons created by two infinitesimally close sources in a far away pulsar 
(or some other source), and travel side by side towards the earth for billions of years. If they interacted at all this 
would clearly have some non-trivial effects on their polarization. This discussion will be fully addressed elsewhere. 


IV. SPIN-SPIN, SPIN-ORBIT, ORBIT-ORBIT INTERACTIONS IN MASSIVE GRAVITY 

Gonsider two conserved sources, = 0, each having the following energy-momentum tensor components [T^ 


Tio 


T^(l) 


T- ■ 


_ t^(2) 

b 5 


7^ _ 7z(0) I rp{2) 

^00 — ^00 + ^00 > 


(29) 
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where the relevant parts read 


^0? = m6{x- Xa), 

'^00 = 7;i^v'^S{x - Xa) - l:J'" v'‘e"''^djS{x - Xa), 


= -mv'^S (x - Xa) + -J^ e^'^^djS {x - Xa ), 
= mv^v^S (x — Xa) + j’‘v^^e^^^''dkS (x — Xa) ■ 


(30) 


Here Xa = Xa(t) is the location of the particle and refers to symmetrization with a 1/2 factor. Note that our 

signature is (—1-++). The relative signs and coefficients are fixed by the requirement that at the first order (namely 
up to 0{v^) and 0{vJ) ) conservation equations are satisfied. Since it could sometimes be a little confusing let us 
explicitly compute two examples. At the first order, 


doT°° + = rriadoS (x - Xa) + rUavldiS {x - Xa) - -Ja (f - fa) 


(31) 


dx 

The last term vanishes obviously, and we have (f — fa) = —VadiS (f — fa) since v = Thus (311 vanishes at 

at 

the first order. Similarly we have 

doT°^+d,T"^ = -mavlv'^dmS {x - Xa)+^JaV'^e>'''^dnd„,5 {x - Xa)+mavlvidi5 {x - Xa)+^Jivlt^'^'^dkdi5 (f- fa) = 0. 

(32) 

In what follows, we will need the following identities 


„ {x'^ - x'^) 

dkv = - -^ = r^ 


r 

_(^fc _ 

dk'r = ----- = -r\ 


1 -(x'^ - x"=) - 


dk- = 


1 {x^ - x'^) 


dk' - = 


r r 

1 1 


dkdn'r = - + r'=f”) , dkO^, - = ^ (5'=" - 


(33) 


where we have assumed, r 0 otherwise one picks up a "Fermi contact" term in the last expression. But this is 
irrelevant for our computation since we are in the long separation regime of two sources. 


A. Gravitomagnetic Effects in General Relativity 


In 71 = 3 + 1, we have 


4 Ut = - 2 KT'oo{d^)~^T^° - 4 Kr'oi(a 2 )-iT°* - 2 KT[^{d'^)-^T^^ + KT'{d‘^)-^T, 
where the trace of the energy-momentum tensor reads 

T = -Too + = -TTlaS (f- fa) + ^Wafa^ (f- fa) “ ^-^a (f “ fa) ■ 

Let us compute the terms separately. The energy density terms read 

Too{d‘^)~^T'°° = miS {x - xi) + ^mivlS (f - fi) - v\T'^^dkS (f - f i) ( 9 ^)“^ 
77126 (^x’ - X2^ + ^m2vi6 {x' - X2^ - ^ J™ ^x' - X2^ 

with the individual terms yielding the following contributions 


(34) 


(35) 


(36) 


c /■-* — 1 X f \ 77117712 

mid (x — xi) (a ) m2d [x' — X2] = - 1 , 

V / 47rr 


( 37 ) 












miS {x — xi) (d^) ^-1x12X26 (^x'— X2^ 


1 77117712^1 

2 47rr 




(38) 


miS (x-xi) (52)-i J 2 ™ (x*' - X 2 ) 


1 777 i(f X V2).J2 

2 47rr^ 




(39) 


1 ^2X1^ 77)21-1 a/-/ -A lmim2i^, 

-miV^d [x - Xi) [d ) 77720 (^a:' - X 2 J = - —t, 

tZ -1/021-1 r / 1 7772(7 X 771 ).Ji 

J{ 77ie (x - Xi) {d^) ^ 1 x 126 yx' - X2j = - - 


4717-2 


t. 


(40) 


(41) 


Observe that we have dropped the higher order term 0{J^v^) to have a consistent expression. The collection of 
these terms yield 


-2KToo(a^)”^T''J‘' = -2 k 


mim2 I ^ vf + vi 
Airr 


1 /777i(f X X 2 ) • J 2 7772(7 X Wl) • Jl 


4777^ \ 


t. (42) 


Similarly the trace-trace interaction reads 
T'{d^)-^T = 


— xni5 {x — xi) +-rriivlS {x — xi) —-j\v\d^^dk6 {x — xi) (d^) ^ 


1x126 ^x' - X 2 ^ -I- ^1x12X26 ^x' - X 2 ^ - ^ v{e^'^'^d'n6 ^x' - X 2 ^ 


which after carrying out the relevant integrals becomes 


KT'{d‘^)-^T = 


777 1 7772 
4717 


1-f 


n ~ '^2 


1 / 7771(7 X V2) ■ J2 7772(7 X Xi) • Jl 


4717 2 


(43) 


(44) 


Observe that T^j{d^) term in (341 yields higher order corrections which we drop, hence the final piece is the 

TQj(9^)“^T°* term which reads 


r',(a2)-iT0* = 


— 777 i 7 i 5 (x — Xi) -|- - Jl d^^dj6 (x — Xi) 


{d^r 


1x12X26 ^x' - X2^ ~ (27' - X2^ 


whose explicit computation gives 

1 x 111 x 12 X 1 •72 11 


-4KT'oi{d^)-^T^^ = -4k 


4„ )-4„» 


(45) 


1 2Ji • J 2 - 3(7 X Jl) • (7 X J 2 ) 


Collecting all these parts above, the potential energy in GR at the desired order becomes 


(46) 


G 

3 _o 3 _2 _ ^' 

G 


UgR = -777 1 7772 

7 

1 + 2^1 + 2^2 - 471 • V2 

fp3 

Jl. J 2 — 3 Jl. 7 J 2.7 


G 

Ji2 


37771(7 X X2) ■ J2 37772(7 X Xi) • Jl 


— 27771(7 X Jl) • J2 -I- 27772(7 X X2) ■ Jl 


(47) 


which include the Newtonian potential energy, plus relativistic corrections such as spin-spin and spin-orbit effects. 
Note that Ji and X 2 are defined with respect to a frame at rest. 
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B. Gravitomagnetic Effects in Massive Gravity 

Analogous computation in massive GR is somewhat more involved but still it is straightforward, we relegate to the 
Appendix the details and quote the final result here. One starts with 

4t/t = -2Kr'oo T - {9^ - - 2kT[^ {9^ - rO (48) 

O 

Again the Tb — T®-! term will not contribute at the order we are working. At the desired order, gravitomagnetic 
potential in mGR is 


UmGR = —Ge 


,mim2 


4 4 ^2 4-2 . - -' 

3 + + 3^®2 - 4ui • V2 


Ge (l + TOgr + m^r^) 


Ji. J 2 — 3 Ji . f J 2 . f 


(1 + nigr + 


Ge-^9r 


(1 + mgv) 


(1 + TOgT + rrigr'^) 

4mi {f X V2) ■ J2 4m2 (f x ui) • Ji 


(49) 


— 2 toi (f X vi) ■ J2 + 2 m 2 (f X U2) ■ Ji 


In the iTLg —>■ 0 limit for not too large distances from (491, one obtains 
UmGR 


-G 


mim2 


4 4 ^2 4 _2 , ^ _ 

3 + - 4vi ■ V2 


G r 

G 

fi2 


Ji. J 2 — 3 Ji. f J 2 . f 


4mi (f X W2) • J 2 4m2 (r x ui) • Ji 


(50) 


— 2 toi (f X vi) ■ J2 + 2m2 (f X 172) • Ji 


So the spin-spin part smoothly reduces to the GR expression in this limit, while a new discontinuity arises (a discrete 
8/9 difference between GR and mGR ) in the O(v^) and 0{vJ) terms. Observe that, if one takes the large r limit 
first, then goes to the nig —> 0 limit, then, as noted in [S], the spin-spin part also gives a distinctly different answer 
from GR. As discussed in the Introduction, a detailed analysis in that work revealed that the discrete difference arises 
for distances > 1.62. 


V. GRAVITOMAGNETIC EFFECTS IN QUADRATIC GRAVITY 

The relevant potential energy is 


nT'id^) V « . 


Uguadt = -'^r^Ad^) + 


+- m})-'- 


KT'(d^-mfj T nr{d^-mfi V 


12 


(51) 


where there are two additional massive modes: a massive spin-2 graviton with ~ ~ A ^ massive spin-0 mode 
with ml = 4 ^( 3 q_|_^) • Massless spin-2 mode of GR is intact and so there will be terms added to (47). Here we define 
Uquad = Ugr + U 2 , where 


U2t = lA, {d- - mj)-^ - m})-^ r°® + (5^ - - ^T' {d^ - mjr T 


, -1 


V -1 , 


. -1 ; 


( 52 ) 
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The ensuing computation is similar so we can simply write the final result 


U2 = 


Gm 1777.2 




G 

^2 


4 4 ^ 

-mi{f X • J 2 - -r772(f X v{) ■ Ji - 2mi(f x hi) • J 2 + 2m2(f x ul) • Ji ) (1 + rm/ 3 ) e 

+ f X V 2 ) ■ J 2 - ^m 2 {r X iTi) • Ji ) (1 + rm^) e”’’"*” 


G 


+ ;:3 i^ + rm^ + r^ml) 


Ji • J 2 — 3 Ji • f J 2 • f 


(^1 + rmp + 

+ rmp + - 


-rm0 


(53) 


Let us compute the r —> 0 limit, where we expect higher curvature terms to play a role, in the potential Uquad- In 
this limit we obtain 

Uquad ^^2 _|_ ^ 2 ^ _|_ constant. (54) 

Observe that all the spin-spin, spin-orbit terms dropped as the quadratic parts cancel the GR parts but a repulsive 
0(v^) term survives. Repulsive nature of quadratic gravity at small separations is expected and in fact, this is the 
reason why the theory is less divergent in the UV regime compared to Einstein’s theory. But it is well-known that 
non-zero /3 gives a spin-2 ghost [T5] . 


VI. CONCLUSIONS 

In our recent work [B], we initiated a study of gravitomagnetic effects in Fierz-Pauli theory (mGR) which is the 
unique viable linear massive gravity that any non-linear version must reduce at the weak field limit (see |14j for related 
issues). In that work we computed the spin-spin interaction between distant sources (such as two galaxy clusters, or 
galaxies) and found that GR and mGR yield different prediction. In the current work, we extended these discussions 
to the case where the sources also have velocities and orbital motion as well as spins. Our computation is based on 
the tree-level graviton exchange diagram and we kept the terms up to and including 0{v^)^ 0{vJ), 0{J^). mGR 
has distinctly different predictions compared to GR in all these orders, extending the well-known vdVZ result in the 
Newtonian potential energy. What is quite interesting is that, especially for the spin-spin interaction part the distinct 
differences arise for separations that satisfy rUg r > 1.62. It is important to realize that for these distances, we are well 
outside the Vanishtein radius m and so the linearized gravity is perfectly legitimate. If one considers the Gompton 
wavelength of the graviton to be at the same order as the Hubble radius of the universe, Ac = ^ — i?„ — 4.4 x 10^®m 
m. then massive gravity shows its most dramatic effects just at the edge of the observable universe. 

Our computations suggest that in addition to the expectation that graviton mass can explain the accelerated 
expansion of the universe, there are further, possibly observable effects of massive gravity on the spin orientations of 
galaxy clusters, or galaxies. It is very early to actually quote data on this but recent observations of unexplained spin 
alignments of 19 pulsars separated for billions of light years m could possibly be due to massive gravity. See also 
the spin alignment observations of galaxies in US], m- 

We have also calculated the gravitomagnetic effects in Einstein -|- quadratic gravity theories, extending the discussion 
of [S] and found that as r —> 0, the potential energy has a repulsive part due to the kinetic energy of sources while 
the Newtonian, spin-spin and spin-orbit interactions vanish. 

Finally, for weak fields and slow velocities, at the lowest order, GR has a nice formulation which resembles Maxwell’s 
theory, where the electric field is replaced by the gravitoelectric field and the magnetic field is replaced by the 
gravitomagnetic field, the electric charge density is the replaced by mass density and the electric current is replaced 
by the mass current. 

For the massive case, the situation becomes more complicated since both in the massive Maxwell’s theory that is 
the Proca theory and massive gravity, Ag and become "physical" hence they must appear in the field equations. 
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VII. APPENDIX: DETAILS OF THE COMPUTATIONS 


Here we give the details of the computation for the massive gravity case. We start with (481 and integrate the 
terms separately. The energy-density interaction part reads 




miS {x — Xi) + {x — xi) — {x — Xi) 


17126 (x' - X 2 ^ + ^7712X^6 (^x' - 0 : 2 ) ~ d '^6 (^x' - £ 2 ) 


= t 


. V^ + V^ \ 1 

mim2 1 H - — ^ -I - — I X 




(271-)- 


(mi J 2 '" + m 2 J{ vle^^'^dk) 


+ - (mivlJ^ + m 24 Ji vle^^'^dk + J( vie^^^Okd'^ 


(27r)2 




(55) 


where the Bessel functions in three dimensions are 


Ki (r nix) = \l 7:,K3 (r nix) = , \/— ( 1 H-| .Ksir nix) = < / — ( 1 H-h 7 -ttt 

2 nix V 2 2 yrrm^V 2 V rnixj = y'r nix V 2 V rnix {rnix)^ 


Now using dz{z (z)) = —z ''K^j^i (z) , we have 


d'k[{rnig) Ki, {r nig)] = —— ^ nig {rnig) '' K^+i{rnig) 

(x^ — x'^) i, 

dk [{rnig)~^Ki, (rnig)] = - ^ - - nig {rnig)~’' K^+i (rnig) 

dkd], [(r nig)~'-'K„ {r nig)] = <5'"” (r nig)~^''^^^ K„+i (r nig) - {r nig)~'' K ^+2 {r nig) 

In order to simplify the computation, we write (|55|) as 


Too{d^-^l} ^ 2 | x (^^g) 2 Ai (rm^) 


+ ^ (mi + m2j( vl^^’^dk 

(27r)5 (2 ^ 




‘ -h ni2vyi vle^^'^dk + J] 

-1 

X \^{rnig) 2 Ki (rnig) 

After carrying out the integrations, one arrives at 

miJ2™i;^e^™X [( rnig) ^ Ki (rmg)l = nigmi{r x V 2 ) ■ J 2 (rnig) ^ Ks (rnig) , 


(56) 


(57) 


(58) 


(59) 


ni2J[v\e'^^dk (rnig) ^ Ki(rnig) =—nigni2{f x vi) ■ Ji (r mg) ^ Ka^rmg). 


(60) 


miUj J™ (rmg) ^ Ki (rmg) = mgmivl{r x V2) ■ J2 (rmg) ^ Ks (rmg) 


(61) 














































12 


1712^2 Jiv\€^^^dk (juig) ^ Ki{rmg) =—mgm 2 V 2 (j x vi) ■ Ji (r nig) ^ Ks (rnig) . 


(62) 


J[ J™ vie^^^dkd'^ (r m,)^ Ki (r m^) 


Collection of these terms gives 


2 

- 



9 

- 


ifl.V2Jl.J2 - Vi.J2V2-Jl 


{rnig) ^ Ki (rnig) 




(27r) = 


— (f X vi).Ji{f X V2).J2 i'f’mg) ^ K 5 (rnig) 


-2 _|_ -2 \ ^ 

miTO 2 ( 1 H-- I (rnig) ^ Ki (rnig) 


+ -mg mi(f X V2) ■ J 2 - ni2{f x Ui) • Ji I (rnig) ^ Ki (rnig) 


Similarly one finds 


(63) 


(64) 




1 m„ 


mini 2 { 1 - ^Vi - ^ 1 ) 2 ] (rnig) ^ Ki{rnig) 


+mg[ - 


6 (27r) 

mi(f X V 2 ) ■ J 2 TO 2 (f X ifi) • Ji 


{rnig) " Ki {rnig) 


(65) 


-/cT', - m2) Vo* = - 


(27r)- 


— niini2Vi ■ V2 {rnig) ^ Ki {rnig) 


-m„ 


mi{f X v{) ■ J 2 ni 2 {r X v‘ 2 ) ■ Ji Ji ■ J 2 


+ "^ 9 — 2 —(^’^g) j{rmg)^ Ki{rmg) 
+ im2(f X Ji) • (f X J 2 ) (rmg)~ K| (rmg) 


( 66 ) 


Collecting all these pieces, one arrives at (49l. 
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